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Abstract 

We have calculated the large-g expansion for the energy cumulants 
and the magnetization cumulants at the phase transition point in the two- 
dimensional g-state Potts model to the 21st or 23rd order in using 
the finite lattice method. The obtained series allow us to give very precise 
estimates of the cumulants for q > 4 on the first order transition point. 
The result confirms us the correctness of the conjecture by Bhattacharya 
et al. on the asymptotic behavior not only of the energy cumulants but 
also the magnetization cumulants for q — > 4+ . 
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1 Introduction 

It is a challenging problem to determine the order of the phase transition in the 
physical system where the order changes from the first order to the second one by 
a small variation of a parameter of the system. It is important to know how the 
quantities that remain finite at the first order transition point diverge or vanish 
when the parameter approaches the point at which the transition changes to the 
second order one. The g-state Potts model[[l], |2| in two dimensions gives a good 
place to investigate the properties of the phase transition around such a point. 
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It exhibits the first order phase transition for q > 4 and the second order one 
for q < 4. Many quantities of this model at the phase transition point j3 = fit 
are known exactly for q > 4, including the latent heat|| and the correlation 
length^, |5|, El, which vanishes and diverges, respectively, in the limit q — > 4 + . 
One the other hand physically important quantities such as the specific heat and 
the magnetic susceptibility at the transition point, which also increase to infinity 
as q — > 4_|_, arc not solved exactly. Bhattacharya et a/.[Q made a stimulating 
conjecture on the asymptotic behavior of the energy cumulants including the 
specific heat at the first order transition point; the relation between the energy 
cumulants and the correlation length in their asymptotic behavior at the first 
order transition point for q — ► 4 + will be the same as their relation in the 
second order phase transition for q = 4 and (3 — > fit, which is well known from 
their critical exponents. If this conjecture is true not only in the g-state Potts 
model but also in other general physical systems with a phase transition whose 
order change from the first order to the second one when some parameter of the 
system is varied, it would give a good guide in investigating such systems. 

Bhattacharya et al. f| also made the large-*/ expansion of the energy cumu- 
lants to order 10 in z = l/y/q and they analyzed the series according to the 
conjecture, giving the estimates of the cumulants at the first order transition 
point for q > 7. The obtained estimates are at least one order of magnitude bet- 
ter than those given by other methods such as the Monte Carlo simulations || 
and the low-(and high-)temperature expansions JTpj, p"i"| , in spite of the fact 
that the maximum size of the diagrams taken into account in their large-g ex- 
pansion is less than the correlation length for q < 15 while the Monte Carlo 
simulations are done carefully on the lattice with its size about 30 times as 
large as the correlation length and also the fact that the maximum size of the 
diagrams in their large-g expansion is as small as the half of the maximum size of 
the diagrams in the low-(and high-) temperature expansions These strongly sug- 
gest the correctness of the conjecture. The large-q series obtained by them are, 
however, not long enough to investigate the behavior of the energy cumulants 
for q closer to 4. 

In this paper we will enlarge the large-g series not only for the energy cumu- 
lants but also the magnetization cumulants at the transition point to order of 21 
or 23 in z using the finite lattice method fl3|, |l4| , |l5| , instead of the standard di- 
agrammatic method used by Bhattacharya et al. The finite lattice method can 
in general give longer series than those generated by the diagrammatic method 
especially in lower spatial dimensions. In the diagrammatic method, one has to 
list up all the relevant diagrams and count the number they appear. In the finite 
lattice method we can skip this job and reduce the main work to the calcula- 
tion of the expansion of the partition function for a series of finite size lattices, 
which can be done using the straightforward site- by-site integration p6| p~7[ with- 
out the diagrammatic technique. This method has been used mainly to generate 
the low- and high-temperature series in statistical systems and the strong cou- 
pling series in lattice gauge theory. One of the purposes of this paper is to 
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demonstrate that this method is also applicable to the series expansion with 
respect to the parameter other than the temperature or the coupling constant. 
Such an approach was successfully done by Bakaev and Kabanovich in genera- 
tion the long series for the ground-state entropy of the g-state antiferromagnetic 
Potts model in two dimensions Q . The long large-g series obtained here by the 
finite lattice method enables us to examine the conjecture by Bhattacharya et 
al. on the asymptotic behavior of the energy and magnetization cumulants for 
q very close to 4. We can also give the estimates of these cumulants for each 
q = 5, 6, 7, • • • that are orders of magnitude more precise than the estimates by 
the previous works and they would serve as a target in investigating this model 
in various contexts, for example in testing the efficiency of a new algorithm of 
the numerical simulation. A brief note on this subject was already presented in 
[ Eg) , where we concentrated only on the specific heat. 

The remainder of this paper is organized as follows. The algorithm of the 
finite lattice method to generate the large-g expansion of the free energy is given 
for the disordered phase in section 2 and for the ordered phase in section 3. The 
relation between the two algorithms is also discussed in section 3. The obtained 
series for the energy and magnetization cumulants at the phase transition points 
are presented in section 4. In section 5 the series are analyzed to investigate 
the asymptotic behavior of the cumulants at q — > 4 + as well as to make the 
estimation of the cumulants at various values of q > 4. In section 6 we conclude 
with a brief summary and a few remarks. 



2 Algorithm in the disordered phase 

The model is defined by the partition function 

Z = J2ex V hj28 SitSj +hJ2S Si ,i), (1) 

\ (iJ) * / 

where the spin variable Si at each site i takes the values 1,2, • • • ,q and (i, j) 
represents the pair of nearest neighbor sites. The phase transition point fit for 
h = is given by exp (/3 t ) — 1 = ^fq. We consider below in this section the free 
energy density in the disordered phase, which is given by 

F d (P)= lim {L x L v )- x hxZ d {L X: L y ) . (2) 

Here the partition function Zd{L x , L y ) in the disordered phase should be cal- 
culated for the L x x L y lattice with the free boundary condition as depicted in 
Fig.l. We first restrict ourselves to the discussion in the case of h = 0. The 
large-g expansion of the partition function in this phase can be given through 
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Figure 1: L x x L y lattice (in this figure L x 
boundary condition for the disordered phase. 



4 and L y — 5) with the free 



the Fortuin-Kasteleyn representation[EO| as 



Z d (L x ,L y ) = q L * L * 



XQA d 



(3) 



with 



4>{X) = (e -1) KX) q< x )- L * L v 

= yb{X) z 2[L*L v -c(X)]-b(.X) 



(4) 



Here is the set consisting of all bonds in the lattice, X is any subset of 
Ad and represents a configuration of bonds, excluding the empty set (the null 
configuration) which gives the leading contribution and is separated in the first 
term of Eq.(g). In Eq.(@) Y = (e? - l)/^/q, z = b(X) is the number of 

bonds in X and c(X) is the number of clusters of sites in X . Two sites connected 
to each other by a bond belong to the same cluster (An isolated single site is 
a cluster). The <f>(X) is called the activity of X. It is useful to note that the 
activity can be rewritten as 



cj>{X) 



X \ b W q -b(X)+l(X) 



Y b(X) z b{X)-2l(X) 



(5) 



where l(X) is the number of independent closed loops formed by the bonds in 
X. 



We define Hd(l x , ly) for each l x x l y lattice (1 < l x < L x , 1 < l y < L y 



as 



H d (l x ,l y ) =ln \Z d (l x J y )/q^y 



(6) 
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where Z d (l x ,l y ) is the partition function for the l x X l y lattice, again with the 
free boundary condition as shown in Fig.l. We next define W d (l x , ly) recursively 
as 

W d (l x , ly) = H d (l x ,l y ) 

- E & + + l)Wd{l' X ,l'y) ■ (7) 

(Oy)^(Wy) 

Here the term with l' x = l x and l' y — l y should be excluded in the summation. 
Then the free energy density defined by Eq.(||) is given by Jl5| 

F d (0)=ln(q) + J2W d (l x ,l y ). (8) 

Now we have to know from what order in z the W d (l x ,l y ) starts. Only 
for this purpose, let us suppose to apply the standard cluster expansion p2f 
based on the diagrams the large-q expansion for the free energy in the disordered 
phase. It should be stressed that in the practical calculation of the finite lattice 
method we need not perform this cluster expansion at all. A polymer is defined 
as a set X of bonds that satisfies the condition that there is no division of X 
into two sets X x and X 2 so that 4>{X) = (f){Xi)<f){X^) . A set X composed of a 
single bond is a polymer. Below in this paragraph we discuss on the polymers 
that consists of two or more bonds. Now we define a cluster of bonds in X 
so that two bonds that are connected by a site belong to the same cluster. 
A polymer consisting of two or more bonds is composed of a single cluster of 
bonds, since a set of bonds composed of more than one disconnected clusters 
of bonds apparently does not satisfy the condition of a polymer. A cluster of 
bonds is a polymer if and only if the cluster of the bonds cannot be divided into 
tow or more clusters of bonds by cutting at any site that is shared by the bonds 
belonging to the cluster. We call such a cluster of bonds as one-site irreducible. 
In fact for a set X r consisting of a one-site reducible cluster of bonds that can 
be divided into clusters X±, X 2 , ■ ■ ■ ,X n of bonds by cutting at a site, we have 
4>{X r ) — <fi(X 1 )<fi(X2) ■ ■ ■ <fi(X n ), which implies that this set of bonds is not a 
polymer. For a cluster of bonds to be one-site irreducible, it is necessary that 
the bonds should form one or more closed loops and that each bond should 
belong to at least one of the loops, which we call closed. An example of a 
closed cluster of bonds is shown in Fig. 2(a). For each closed cluster of bonds, 
the loop of the bonds that forms the perimeter of the whole cluster is called 
the exterior boundary. The exterior boundary for the closed cluster of bonds in 
Fig. 2(a) is indicated in Fig. 2(b). An interior boundary is a loop inside which 
all bonds are removed and we call the part of the cluster surrounded by the 
interior boundary as a hole. The closed cluster of bonds in Fig. 2(a) has a hole 
with the interior boundary, as indicated in Fig. 2(c). For a cluster of bonds to be 
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Figure 2: (a) Example of a closed cluster of bonds with a hole, (b) The same 
as (a) with the exterior boundary highlighted by thick lines, (c) The same as 
(a) with the interior boundary highlighted by thick lines. 




Figure 3: Example of a polymer whose exterior boundary has a crossing. 



one-site irreducible, not only should it be closed but also its exterior boundary 
has no crossing. If the exterior boundary has a crossing, as in Fig. 3, the cluster 
of bonds is one-site reducible and is not a polymer. Thus a polymer is the closed 
cluster of bonds with its exterior boundary having no crossing. We list as an 
illustration in Fig. 4 all the polymers that can be embedded within 3x3 lattice 
with their activities. The polymers are skipped that can be obtained by rotating 
and/or reflecting those in Fig. 4. For later convenience we add some definitions 
concerning a polymer. A concave corner in the boundary of a polymer is the part 
of the boundary composed of two neighboring bonds that are bent perpendicular 
to each other toward the interior of the polymer. If two neighboring corners on 
the boundary are both concave, we call the two corners together with the bonds 
in a straight line connecting the two corners as a concave part of the boundary. 
Examples of the concave parts are given in Fig. 5. 

According to the theorem of the standard cluster expansion |2l[ |22j, the free 
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Figure 4: All the polymers that can be embedded within 3x3 lattice with 
their activities. 

energy for the l x x l y lattice is given by 
\n.Z d (l x ,l y ) = lj y laq 

+ c(X u X 2 ,---,X n )4>(X 1 )ct>(X 2 )---4,(X n ), (9) 

{Xi ,Xi ,-■ ■ ,X n -.connected} 

where the summation should be taken for every cluster of connected polymers 
{Xi, Xi, ■ ■ ■ , X n } defined by the condition that there is no division of the cluster 
of polymers {X\, X 2 , • • • , X n } into two clusters of polymers {X[, X 2 , • • • , X' m } 
axid{X'^,X! 2 ',---,X^_ m } (to < n) so that <f)(X 1 \JX 2 \J ■■■ UX n ) = c/){X[UX 2 U 
■ • • U X' m )(j)(X'{ Ulj'U ■ • • U X'n_ m ). (The 'cluster' of polymers here should be 
distinguished from the 'cluster' of sites in the Fortuin-Kasteleyn representation 
or the 'cluster' of bonds defined in the explanation of the polymers above.) The 
c(Xi, X 2 , • • • , X n ) in Eq.(^|) is the numerical factor that depends on the way how 
the clusters of polymers are connected, the details of which are not necessary 
here. 

When a cluster of connected polymers can be embedded into the l x x l y 
lattice but cannot be embedded into any of its rectangular sub-lattices, we call 
it shortly to fit into the l x x l y lattice. Then we can prove the theorem that 
the Taylor expansion of the Wd(l x , ly) with respect to z includes the contribution 
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Figure 5: Example of a polymer with two concave parts highlighted by thick 
lines. 

from all the clusters of polymers in the standard cluster expansion that fit into 
the l x x l v lattice. 

Now, according to this theorem, we can know the order of Wd(l Xl l y ) hi z 
by examining the clusters of connected polymers that ht into the l x x l y lattice. 
We first note that Wd(l,l) = since the lxl lattice consists of a single 
site with no bond, so that no polymer can be embedded into this lattice, and 
Wd(2, 1) = Wd(l,2) = 0(z) since the 2 x 1 or 1 x 2 lattice consists of two 
sites connected by one bond, so that a polymer composed of a single bond can 
be embedded into this lattice. We also note that Wa(l x , 1) = Wd(l, l y ) = for 
l x > 3 or l y > 3, because the corresponding lattice consists of the sites connected 
by bonds all on a straight line and the only clusters of polymers that fit into 
the lattice are those with each polymer composed only of a single bond and 
these polymers are not connected according to the definition of the connected 
polymers. 

We next consider the order of Wd{l xi l y ) for l x > 2 and l y > 2. We first 
consider a cluster consisting of a single polymer that fits into the l x x l y lattice. If 
the polymer has a concave part in its exterior boundary we can add some bonds 
to make the part flat, obtaining a new polymer whose activity is changed by a 
factor of Y n z~ 1 with n the number of added bonds, so its order in z is lowered. 
This implies that a polymer that has concave part in the exterior boundary does 
not contribute to the lowest order term of Wd(l x ,l y ) in z. Let us next consider 
the case when the polymer has some holes. The interior boundary surrounding 
each hole has at least one concave part. (We take the hole as the exterior of 
the polymer when we say 'concave' here.) We can add some bonds to make the 
concave part flat, again obtaining a new polymer with the order of its activity 
in z lowered by a factor of z^ 1 . This implies that a polymer that has a hole does 
not contribute to the lowest order term of Wd(l x , l y ) in z. Thus the remaining 
polymers that have neither a concave part on the exterior boundary nor a hole in 
the interior give the contribution to the lowest order term of Wd(l x ,l y ), examples 
of which are shown in Fig. 6. They can be converted into each other by adding or 
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Figure 6: Examples of the cluster consisting of a single polymer that contribute 
to the lowest order term of the Wd{l x , ly) with l x = 4 and l y = 5. 

removing pairs of bonds to their concave corners or from their convex corners, 
respectively. When each pair of bonds is added or removed the activity of the 
polymer is changed by Y 2 with its order in z unchanged. Hence we see that 
their orders are all the same; that is, the order of z lx+l y~ 2 . Next we consider a 
cluster of connected tow or more polymers that fits into the lattice. It is enough 
to take into account the clusters of polymers that do not have either a concave 
part on the boundary or a hole in the interior, since we are going to know the 
clusters of connected polymers that contributes to the lowest order term in z. 
Among the clusters of connected two or three polymers that fit into the lattice, 
those give the lowest order contribution in which the polymers touch each other 
only through their exterior boundary and each neighboring pair of polymers 
share only one bond, like Fig. 7(a), and this type of clusters have the order of 
z l x +ly+n-3 £ or j-^g cnis t cr f n polymers (n — 2,3), which is higher than the 
lowest order contribution from the clusters composed of a single polymer. Wc 
note that a cluster of two or three polymers in which each neighboring pair of 
polymers share only one site, an example of which is given in Fig. 7(b), is not 
connected according to the definition of the connected polymers. Among the 
cluster of connected more than four or more polymers that fit into the lattice, 
those give the lowest order contribution in which four of the polymers surround 
a plaquette and each neighboring pair of polymers share only one site, like Fig. 8, 
and this type of clusters have the order of z l ^+ l v+ n ~ i (n > 4) for the cluster 
of n polymers, which is also higher than the lowest order contribution from the 
clusters composed of a single polymer. In summary the lowest order term of 
Wd{l x >ly) has the order of z l:c+l y~ 2 for l x > 2 and l y > 2, which comes from 
the clusters of a single polymer without a hole or a concave part fitting into 
l x x l y lattice. Thus in order to obtain the series to order z N for the free energy 
density, we have only to take into account in Eq.(0) all the rectangular lattices 
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(a) (h) 

Figure 7: (a) Example of the cluster consisting of connected two polymers that 
share only one bond, (b) Example of the cluster consisting of disconnected two 
polymers that share only one site. 

that satisfy l x + l y — 2 < N. 

We summarize the procedure of the finite lattice method to generate the 
large-g expansion series of the free energy density in the disordered phase. 

1) Calculate the partition function Zd(l x ,L y ) defined by Eq.(||) as a power 
series of z to order z N (and of Y in full order) for each l x x l y lattice with 
l x + l y — 2 < N. In these calculations use the transfer matrix method and inte- 
grate bond by bond|ll| [l7|], which reduces the necessary CPU time and storage. 

2) Calculate Hd{l x J y ) defined by Eq.(^|) as a power series of z by Taylor ex- 
panding \nZd(l x , l y ) for each l x x l y lattice. 

3) Calculate Wd{l Xl l y ) defined recursively by Eq.(Q) for each l x x L lattice. 

4) Sum up all the Wd(l x , ly)'s with l x + l y — 2 < N according to Eq.(||) to obtain 
the expansion series for the free energy density. 

Now we discuss on the free energy in the case of h ^ 0. The Boltzman factor 
in Eq. ([!]) that involves the effect of the magnetic field h can be divided into two 
parts as 

exp (hS Sii i) = 1 + [i5 Si ,i (10) 

with n = e h — 1 . The partition function Zd(l x ,l y ) is given by the same expression 
as Eq. (||) with X now representing the configuration of both bonds and frozen 
sites. These frozen sites in X are the sites whose spin variables are fixed to 1 
corresponding to the second term of Eq. (|l0|) . The activity </>(X) is now given 
by the right hand side of Eq.(^|) multiplied by the factor /j^W with s(X) repre- 
senting the number of the frozen sites in X and by the factor q^^ x ^> = z 2 ^ x ^ 
with f{X) representing the number of the clusters of sites that include at least 
one frozen site in X . Hd(l x , l y ) and Wd(l x , l y ) are defined by the same equation 
as Eq.(^|) and (Q) and the free energy density in the thermodynamic limit is 
given by Eq.(^[). In the standard cluster expansion the polymer whose activity 
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Figure 8: Example of the cluster consisting of connected four polymers that 
surround a plaquette and each neighboring pair of polymers share only one site. 

is proportional to /i" is given by the cluster of bonds and n frozen sites that 
cannot be divided into two or more clusters by cutting at a site so that one of 
the clusters involves no frozen site. There is only one polymer whose activity 
is proportional to /i, which is composed of a single frozen site with no bond. 
Therefore, if we expand Wd{l x , ly) m A* as 

oo 

W d (l x ,ly)=J2 W d n) (lM /»" , 
n=0 

we know that W^\l x , l y ) = z 2 for l x = l y = 1 and W^\l x , l y ) = for the other 
pairs of l x and L. So the coefficient of the term proportional to \x in the free 
energy density is just z 2 . This will lead to the result that the magnetization 
in the disordered phase is zero. We can see for n > 2 that wf} (l x , l y ) has 
the order of z( iar +*») for n < l x l y and has the order of z k ( l *+ l y) or higher with 
k = [n/ (l x l y )] for n > l x l y . Here [p] indicates the integer that is less than p. The 
cluster of polymers that contributes to its lowest order term is composed of a 
single polymer that has n frozen sites and fits into l x x l y lattice for n < l x l y and 
it is composed of k polymers (k > 1) each of which has l x l y frozen sites and fits 
into l x x l y lattice and a polymer with n — kl x l y frozen sites for n > l x l y . Thus 
the order of W^ n \l x ,l y ) for n > 2 is z ll:+ly or higher. So, in order to generate 
the large-g expansion of the term in the free energy density that is proportional 
to fi n with n > 2, it is enough to take into account the finite-size lattices that 
satisfy l x + l y < N in Eq.(§). 
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Figure 9: L x x L y lattice (in this figure L x = 4 and L y = 5) with the fixed 
boundary condition for the ordered phase. The spins of the sites on the bound- 
ary (open circles) are fixed to a single value among 1, 2, • • • , q. 



3 Algorithm in the ordered phase 

Next we consider the large-q expansion of the free energy density F ((3) in the 
ordered phase, which is given by the same equation as Eq.(Q) with the suffix 
d replaced by o where the partition function Z (L X , L y ) in the ordered phase 
should be calculated for the L x x L y lattice with the fixed boundary condition 
as depicted in Fig. 9. The values of the spins on the boundary are all fixed to a 
single value among 1, 2, • • • , q. Here we take this value as 1. 

We first restrict ourselves to the discussion in the case of ft. = again. By 
the duality the ordered and disordered free energy densities are related as 

FM =F G8) + 21n((^-l)A/?) 

for (e? - l)(e? - 1) = q . (11) 

(n) 

One method to obtain the large-q series for the energy cumulant F at p = pt 
in the ordered phase is to use the relations 

2(1 + 2) , 

-3z[FP +FW} + (z-z 2 )[FV+ F p ]} 
-6z[fP-F<;V}-(z-6 Z *)[fP-FM], 
-\0z[Ff + fW] - (5z - 60z 3 )[^ 2) + F&] 
+(z -z 2 - 24z 3 + 24z 4 )[i^ 1) + fW] , 
-15z[Fj 5) - F Q ( 5 )] - (15* - 300^ 3 ) [Fj 3) - F&] 



(i) 

i 

(2) 
i 

(3) 



x o 
7?(2 

O 

^(3 



F (4) _ p (i 



F 



(5) 
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-{z - $Qz A + 360z b )[F^ 1] - F^ L >] 



n(i)i 



(12) 



which are derived directly from the duality relation (p"l|). These were given in 
Ref.j8| except for the last one. 

In spite of that, we give below the algorithm to generate the series for the free 
energy density in the ordered phase by the finite lattice method without relying 
on the duality relation. This algorithm will prove important when we extend our 
work to generate the magnetization cumulants, since they can be obtained by the 
derivative of the free energy density with respect to the external magnetic field 
and we have no longer duality relation for the free energy density in the presence 
of the external magnetic field. The large-g expansion of the partition function 
in this phase can be given also through the Fortuin-Kasteleyn representation as 



Z {L X , L y ) 



= (e? 



X 



(13) 



with 



v&) = (^-ir bp v*> 

= yb(X) z b(X)-2c(X) _ 



(14) 



Here X is a configuration of the removed bonds and we also use X to represent 
the set of removed bonds. In the leading configuration for the fixed boundary 
condition all the pairs of nearest neighbor sites are connected by the bonds 
and there is no removed bonds (X is the empty set), and its contribution is 
separated as the first term of Eq.@. In Eq.([j~4|) Y = ^/(e' 3 - 1) and b(X) 
is the number of bonds in X and c(X) is the number of clusters of sites in the 
configuration of bonds with X removed without taking into account the clusters 
including the sites on the boundary. We define H a (l x , l v ) for each l x x l y lattice 
(1 < l x < L x , 1 < l y < Ly) as 



H (l x ,l y ) = In Z (l x , l y )/ ( 



O0 



1 



(15) 



where Z (l x ,l y ) is the partition function for the l x x l y lattice, again with the 
fixed boundary condition as shown in Fig. 9. We next define W (l x ,l y ) recur- 
sively in the same way as in Eq.(Q). Then the free energy density in the ordered 
phase is given by 



F (p) ^2\n(e l3 -l) + J2 ly) 



(16) 



It is useful to compare now the large-g expansions by the finite lattice method 
for the disordered phase presented in the previous section and for the ordered 
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Figure 10: (a) Example of the configuration X of bonds (solid lines) for the 
free boundary condition, (b) The configuration X of removed bonds (crosses) 
for the fixed boundary condition with X ~ X. 



phase presented here. We first note that the lattice with the free boundary 
condition in Fig.l and the lattice with the fixed boundary condition in Fig. 9 
are dual to each other in the sense that each site of one lattice is located at the 
center of a plaquette of the other lattice, when the two lattices have the same 
size. We denote X ~ X as the relation between the configuration X of the 
bonds on the lattice with the free boundary condition and the configuration X 
of the removed bonds on the lattice with the fixed boundary condition, if the 
location of each bond (identified with its center) in X coincides with the location 
of each corresponding removed bond (also identified with its center) in X (see 
Fig. 10(a) and (b)). Then, the number l(X) of independent closed loops for each 
configuration X in Eq.(^) is equal to the number c(X) of clusters of sites in the 
configuration with X removed in ( |l4| ) for X ~ X (see Fig. 10(a) and (b)), so 
that (j>{X) = ip(X) when Y — Y. We therefore know from Eq.(|), (g), (§J) and 
( p"5| ) that Hd(l x , l y ; Y, z) — H (l x , l y ; Y, z) when Y = Y . hence the same is true 
for Wd(l x , ly\ y, z ) and W (l x ,l y ; Y, z). Thus W (l x , l y ) has the same order in z 
as Wd(l x ,l y ), that is, the order of z lx+ly ~ 2 . So in order to generate the series 
to order z N we have only to take into account the rectangular lattices with 



l x + l y — 2 < N in Eq.(|lg|). This conclusion can also be reached by repeating 
the discussion based on the standard cluster expansion using the diagrams. The 
discussion goes completely in parallel with those given for the disordered phase 
in the previous section if the set of bonds X representing the configuration of 
bonds for the disordered phase is replaced by the set of removed bonds X for 
the ordered phase. 

The relation of Wd(l x , l y ; Y, z) = W (l x , l y ; Y, z) for Y = Y also implies that 
the free energy densities i*d(/3) given by Eq.(||) and F o {0) given by Eq.([l6|) 
satisfy the duality relation ([ill), because the condition that Y — Y is equivalent 
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to the condition for (3 and (3 in the duality relation in Eq.(p"l|). 

Now we proceed to the discussion of the free energy in the case of h ^ 0. 
The partition function Z (l x , ly) is now given by the same expression as Eq.([l3|) 
with X representing the configuration of both removed bonds and frozen sites. 
The activity 4>{X) here is given by the right hand side of Eq.(p^) multiplied by 
the factor /i s W with s(X) representing now the number of the frozen sites in 
X and by the factor q~f( x ^ — z 2 ^ x ^ with f(X) representing the number of the 
clusters of sites that include at least one frozen site in X without taking into 
account the clusters including the sites on the boundary. In the standard cluster 
expansion the polymers whose activity is proportional to u n (n > 1) are given 
by the same polymers that appear in the cluster expansion in the case of h = Q 
discussed above except that their n sites are frozen. These frozen sites should 
belong to the cluster of sites that does not include the site on the boundary. 
By freezing n sites the order of its activity in z increases by In. There is one 
exceptional polymer. It is composed only of a single frozen site with no removed 
bonds and its activity is just a. If we expand W {l x , l y ) in \x as 

oc 

W {l x ,l y ) = Y,Wi n \l x ,ly) M", 
71=0 

we know that Wo^ (l x , L) = for l x — 1 or l y — 1 since such a lattice does 

not have a site except on the boundary. We also know that Wj (l x , ly) = 0(1) 
for l x = l y = 2. The cluster of polymers that contributes to its lowest order 
term is composed of n exceptional polymers each of which are composed only 
of a single frozen site. We can see that wi n \l x , l y ) for the other pairs of l x and 
l y has the order of z (^+'«- 2 ). The cluster of polymers that contributes to its 
lowest order term is composed of a polymer that has no frozen site and fits into 
the lattice with the order of z(^+ ; »- 2 ) and n exceptional polymers each of which 
are composed only of a single frozen site. So, in order to generate the large-g 
expansion of the term in the free energy density that is proportional to /i™, it is 
enough to take into account the finite-size lattices that satisfy l x + l y — 2 < N 
in Eq.dJ). 

4 Series 

The algorithms in the previous two sections give the large-g expansion of the free 
energy density at the arbitrary inverse temperature (3. The energy cumulants 
can be obtained by taking its derivative with respect to (3. If we are interested 
in the energy cumulants only at the phase transition point, as is the case here, 
we can set Y = 1 + y and Y — 1 + y (with y = and y — at (3 = /3 t )|2j|, then 
in order to obtain the nth cumulant we have only to keep the expansion with 
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respect to y to order y n as 



F, 



(n) 



P=Pt 



E 



(n,m) 



(17) 



using ^ = (1 + y ■ 



z)4- and 

) dy 



p( n ) 



d n 
1/3™ 



Fa(P) 



P=0t 



E 



(n,m) m 



(18) 



USing ± = [l + y+(l + y)2 2 ]^. 

Using this procedure we have calculated the series to order N = 23 in z for 
the nth cumulants with n = 0, 1, 2 and to order ./V = 21 for the cumulants with 
n = 3, ■ • • , 6. The coefficients of the series are listed in Table 1 and 2 for the 
disordered phase and in Table 3 and 4 for the ordered phase. We have applied 
the finite lattice method both in the ordered and disordered phases without 
relying on the duality relations ( [l2|) for the cumulants. We have checked that 
all of Wd{l x , l y ) and W (l x ,l y ) with l x + l y — 2 < N have the correct order in 
z as described in the previous two sections. The obtained series for the zeroth 
and first cumulants (i.e., the free energy and the internal energy) agree with 

In) 

the expansion of the exactly known expressions. The coefficients for F G agree 
with those by Bhattacharya et al. to order 10. We have also checked that the 
series satisfy the duality relations ( |l2| ) for the cumulants. 

The nth magnetization cumulant of the free energy at (3 = j3t and h = is 
defined by 



1 



d" 



K q-l dh r < 



F 



d , o 



l 



/i=0,/3=A 



9-1 



E 6 £ m) * r ' 



(19) 



It can be calculated from the expansion series of the free energy in n to order 
u" using jfc = (1 + We have calculated the series for the nth cumulants 

with n = 1, 2, 3 to order N = 21 in z. The coefficients of the series are listed in 
Table 5 and 6 for the disordered and ordered phases, respectively. The series for 
the first cumulant (i.e., the magnetization) in the disordered phase is vanishing 
and not listed here. The series for the first cumulant in the ordered phase agrees 
with the expansion of its exactly known expression. 
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Table 1: The large-g expansion coefficients a^™ for the nth energy cumulants 
of the free energy at = fit in the disordered phase (n = 1-4). 



m 


(l.m) 
a d ' 


VI, m) 
a d ' 


(3.m) 

a d ' 


(4,m) 
d 



1 


2 


2 


2 


2 


2 


4 


14 


58 


242 


3 


-2 


26 


338 


3026 


4 


2 


118 


2474 


37186 


5 


-6 


250 


11490 


291946 


6 


4 


894 


55506 


2068914 


7 


-16 


1936 


219404 


12144136 


8 


6 


6160 


876384 


65992348 


9 


-38 


13538 


3127818 


322997186 


10 




39774 


11161914 


1494569130 


11 


-76 


88360 


37125056 


6470376688 


12 


-56 


245188 


122847224 


26851465312 


13 


-96 


547468 


387657840 


106378391164 


14 


-348 


1457976 


1214792028 


407805474096 


15 


156 


3264012 


3679570740 


1510368315660 


16 


-1634 


8410284 


11061186460 


5448643138308 


17 


1946 


18868858 


32420680250 


19136730128506 


18 


-6852 


47391870 


94306690802 


65779858359354 


19 


10744 


106180532 


269004235012 


221297704673276 


20 


-27004 


261607968 


761823806996 


731233034945892 


21 


48492 


586199668 


2124082332388 


2373735516063508 


22 


-102226 


1415498882 






23 


199386 


3174300542 
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Table 2: The large-*? expnsion coefficients a^™ for the nth enrgy cumulants 
of the free energy at (3 = fit in the disordered phase (n = 5,6). 



m 


(5,m) 


(6,m) 
a d 



1 


2 


2 


2 


994 


4034 


3 


24338 


186626 


4 


488162 


5995138 


5 


5985042 


110362330 


6 


60977034 


1583071794 


7 


503300684 


17852715736 


8 


3668770296 


171621163420 


9 


23690762442 


1432076136578 


10 


140396232786 


10720828033434 


11 


768231161504 


72985931413600 


12 


3949989951344 


459566465021728 


13 


19184686377744 


2701568006995708 


14 


88891085678652 


14974796851686096 


15 


394525773133188 


78777820375010652 


16 


1687573606453276 


395842580467193364 


17 


6978567585546266 


1908650742049266298 


18 


28012258546890218 


8869906520648194890 


19 


109405587397512484 


39863282911219604012 


20 


416940286121010524 


173802434948834045508 


21 


1553287387425965092 


737013278292605504308 
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Table 3: The large-g expansion coefficients Oo ' m for the nth energy cumulants 
of the free energy at = fit in the ordered phase (n = 1-4) . 



m 


(l,m) 


(2,m) 

a> 


(3,m) 

ah 


(4,m) 



1 


2 








2 


-4 


16 


-64 


256 


3 


2 


34 


-430 


3778 


4 


-2 


114 


-2654 


41778 


5 


6 


254 


-12186 


322670 


6 


-4 


882 


-57018 


2210982 


7 


16 


1944 


-224732 


12819264 


8 


-6 


6128 


-888024 


68657204 


9 


38 


13550 


-3164682 


333583598 


10 




39698 


-11243178 


1532324246 


11 


76 


88360 


-37363472 


6604807168 


12 


56 


245036 


-123377384 


27298396784 


13 


96 


547356 


-389128512 


107855738700 


14 


348 


1457784 


-1218076500 


412466192928 


15 


-156 


3263316 


-3688318020 


1524965526804 


16 


1634 


8410596 


-11080768444 


5492850472332 


17 


-1946 


18865590 


-32471142890 


19269581858838 


18 


6852 


47395762 


-94419894146 


66169209300214 


19 


-10744 


106166828 


-269288597908 


222430064188868 


20 


27004 


261629456 


-762460849076 


734462791941548 


21 


-48492 


586145660 


-2125652044660 


2382881224028156 


22 


102228 


1415594740 






23 


-199212 


3174081000 
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Table 4: The large-g expansion coefficients a„°m for the nth energy cumulants 
of the free energy at (3 = (3 t in the ordered phase (n = 5,6). 



m 


(5,m) 
Ho 


(6,m) 
Ho 



1 

2 


-1024 


4096 


3 


-29518 


218914 


4 


-556382 


6813714 


5 


-6773994 


126069374 


6 


-67122114 


1774583142 


7 


-546094604 


19774354944 


8 


-3918393456 


187361651588 


9 


-25037212842 


1545876302510 


10 


-146961943266 


11451708807878 


11 


-798499755344 


77296110810160 


12 


-4080741048224 


483066658347536 


13 


-19726182681504 


2822025943834956 


14 


-91033421848212 


15558449192797824 


15 


-402728474968788 


81476650613568516 


16 


-1717926911825116 


407801271671111676 


17 


-7087982430811466 


1959732852092226390 


18 


-28396321463796218 


9080903533904186662 


19 


-110725077538793524 


40709407799692233428 


20 


-421377562607301884 


177104383529723970956 


21 


-1567944342623768692 


749588012269824549980 
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Tabic 5: The large-g expansion coefficients b ™' m for the nth magnetization 
cumulants of the free energy at (5 = (3 t in the disordered phase (n = 2-3) 



m 


,(2,m) 


, (3,m) 
°d 



1 






2 


1 


1 


3 


4 


12 


4 


21 


120 


5 


72 


804 


6 


289 


4943 


7 


920 


25764 


8 


3197 


126310 


9 


9632 


564780 


10 


30589 


2412549 


11 


88412 


9714144 


12 


264037 


37758464 


13 


738508 


140800740 


14 


2109461 


510552463 


15 


5744540 


1794840960 


16 


15861797 


6167896002 


17 


42249692 


20689927908 


18 


113575441 


68107803125 


19 


296922168 


219877596348 


20 


781015901 


698663386396 


21 


2009401452 


2184569257680 
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Tabic 6: The large-g expansion coefficients b " for the nth magnetization 
cumulants of the free energy at /? = (3t in the ordered phase (n = 1-3) 



n 




A2,m) 
<Jo 


1,(3, m) 



1 


1 






2 


-1 


1 


-1 


3 




4 


-12 


4 


-3 


23 


-122 


5 




80 


-828 


6 


-9 


327 


-5159 


7 




1048 


-27156 


8 


-27 


3649 


-134398 


9 




11020 


-605496 


10 


-82 


34902 


-2604438 


11 




100840 


-10546968 


12 


-254 


299901 


-41207005 


13 




837240 


-154329576 


14 


-804 


2380574 


-561777935 


15 




6465720 


-1981489368 


16 


-2596 


17772083 


-6829373776 


17 




47195596 


-22967338956 


18 


-8520 


126320724 


-75775684452 


19 




329203956 


-245115502776 


20 


-28311 


862396671 


-780219055115 


21 




2211793232 


-2443319594304 
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5 Analysis 



In the first part of this section we analyze the large-g series for the energy 
cumulants presented in the previous section. The series of the magnetization 
cumulants are analyzed in the rest part of this section. The asymptotic behavior 
of the energy cumulants at j3 = fit for q — * 4 + was conjectured by Bhattacharya 
et as follows. It is well known that there is an asymptotic relation between 
the energy cumulants and the correlation length £ for j3 — > fit at q = 4 as 

ff \ ~ A T( £-^ e n/2 -\ (20) 

with the constant A that is independent of n. This relation comes from the fact 
that the critical exponents of the correlation length and the specific heat in the 
second order phase transition are v — 2/3 and a = 2/3, respectively, for q = 4. 
Their conjecture claims that the relation ( |20| ) is also held as the asymptotic 
relation between the energy cumulants and the correlation length at the first 
order phase transition point for gr — > 4 + . Then, from the known asymptotic 
behavior of the correlation length at the first order phase transition point for 
q^A + asfj 

f-^x, (21) 

where x = exp with 9 defined by 2 cosh 6* = ^fq (9 ~ *Jq — 4/2 for q — > 4 + ), 

the conjecture implies that the asymptotic form of the energy cumulants in the 
limit of q — > 4 + is given by 

F W j ^ aB n-2 T (n-li) x 3n/2-2 (22) 

We notice that the constants a and i? should be independent of n and from the 
duality relations (|l2] ) that they should be common to the ordered and disordered 
phases. 

To analyze the large-g series for these quantities, Bhattacharya et al. used a 
sophisticated method. Let f(z) be a function that has the asymptotic behavior 
of the right hand side of Eq.(|22|). Then Olnf(z) is a regular function of z 
except for the possible terms proportional to the inverse powers of x, which 
vanishes rapidly when 9 approaches zero. Although 9 itself cannot be expanded 
as the power series of z since 9 diverge at z = 0, the quantity 1 — e~ e = 
1 — 2z/(vl — 4z 2 + 1) behaves like 9 for 9 — > and can be expanded as a 
power series of z. Then one may expect that (1 — e~ e ) ln/(z) will be a smooth 
function of z and its Pade approximant will give a convergent result for q — > 4+. 
In fact using the 10th order series Bhattacharya et al. found that the Pade 
approximants of F^ 2 ' give convergent results for q J> 7. The obtained results 
were more precise than the results of the Monte Carlo simulations H and the 
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low-(or high-)temperature series|l(], O, [l2| as mentioned in the introduction 
and exhibits rather flat behavior of F^'/x around 7 < q < fO, suggesting the 
correctness of the asymptotic behavior (^) for n = 2 with a ~ 0.76. If we 
use our longer 23rd order series of to repeat these Pade analyses, we find 
that the convergent region of q is extended down to q ~ 5, with the accuracy 
of a few percent at q — 5 and the data of jx is rather flat for 5 < q < 10, 
which gives a stronger support to the conjecture that F^ is approaching its 
asymptotic form J2^). 

We, however, adopt here another method of the analysis that is simpler and 
more efficient. The latent heat £ has the exact asymptotic form as (7| 

so, if has the asymptotic form of Eq.(|2|), the product F^C P is a smooth 
function of z for q — > 4 + when p — 3n — 4, while it behaves like some positive 
or negative powers of x and has an essential singularity with respect to z at 
q = 4 when p 7^ 3n — 4. Thus if we try the Pade approximation of F^C P it will 
lead to a convergent result for p — 3n — 4, while the approximants will scatter 
when p is off this value. In the Pade approximation of F^C P the series for 
the ordered and disordered phases should be treated in a slightly different way, 
corresponding to the difference of the lowest order of the series; the series for 
the ordered phase starts from the order of z 2 , while the series for the disordered 
phase starts from the order of z. Thus we should solve the following equation: 

p(n) c 3n-i = z rp M ( z )/Q L ( z ) + Q^M+L+r+lj ^ ^ 

where r is 2 for the ordered phase and I for the disordered phase, and Pm{z) 
and Ql{z) are polynomials of order M and L respectively with M + L + r = N 
where N is the truncated order of the series. Fig. 11 is the result of this analysis 

(2) 

of F^ C p for q = 4. We see very clearly that the Pade approximants have 
the best convergence at p — 2.00(1) and the convergence becomes bad rapidly 
when p leaves this value. All the other approximants for F^j C p have similar 
behaviors with respect to p with the best convergence at p = in — 4, although 
we do not show them in figures. These results are enough to convince us that 
F^ n J has the asymptotic behavior that is proportional to x 3n / 2 ~ 2 . 

Then we calculate the amplitude of F^ for each q by taking the Pade ap- 
proximants of p( n )£3n-4 anc j mu itiplying £ 4-3 ™ that is calculated from its exact 
expression. Usually a diagonal approximant (that is, M = L) or its neighbor- 
hood produces good estimates. We, however, adopt here all the approximants 
with M > 8 and L > 8, since we have found in this case that each of these 
approximants is scattered almost randomly around their average value and we 
have found no reason to give a special position only to the diagonal approxi- 
mant. Some denominators of the approximants have zero at z < 1/2 (q > 4) 
and we have excluded such approximants. The results for n — 2, • • • , 6 in the 
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(2) 

Figure 11: Pade approximants or 'LP at q = 4 with M > 9 and L > 9 
plotted versus p. 

ordered and disordered phases are listed in Table 7-10 for integer values of 
q > 5, with the second energy cumulant transformed to the specific heat by 
multiplying 01 . In Figs. 12-16, we also plot \F^\/x~~ 2 versus 9 for n =2-6, 
respectively. In these figures the dotted and dashed lines indicate the errors for 
the disordered and ordered phases, respectively. In these tables and figures the 
given errors represent the magnitude of the fluctuation in the data of the Pade 
approximants. 

In any case, the estimates have surprisingly high accuracy. Our data are 
orders of magnitude more precise than the result of the Monte Carlo simulation 
by Janke and Kapplcr[^| for q = 10, 15, 20 (also listed in Table 7-10), which was 
performed carefully in large enough lattices checking the finite size effects, and 
they are consistent with each other except for small differences in the data of 
the 5th and 6th cumulants for the ordered phase. On the other hand, the results 
of the low-temperature series [[[2) (also listed in Table 7 and 8) are inconsistent 
with these. This would be because we do not have the correct knowledge (or 
conjecture) of the possible singularity structure of the cumulants at some f3 < fit 
in extrapolating the low-temperature series toward (3 = t . We should empha- 
size that the value of the specific heat has an accuracy of about 0.1 percent 
at q = 5 where the correlation length is as large as 2500 lattice spacingQ and 
there is no previous data to be compared with, and of a few percent even at 
the limit of q = 4 where the correlation length diverges. The accuracy is better 
by a factor of more than 10 when compared with the result of using the Pade 
approximation for (1 — e~ e ) \nF^ n \ The reason for this is the following. Let 
us suppose that the Pade approximants give estimates of F( n )£ 3n ~ 4 = c (i ± S) 
and (1 - e- )\nF^ = c'(l ± 8'), then the estimate of F^ from the former 
would also have an ambiguity factor of 1 ± S, while the estimate of F^ from 
the latter would have an ambiguity factor of e ±c 5 ^ e , which grows up for the 
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Table 7: Estimates of the specific heat Cd for the disordered phase in comparison 
with the low-temperature series and the Monte Carlo simulations. 



q 


large-g 


low-temp. 


Monte Carlo 


5 


2889.1(41) 






G 


205.930(55) 






7 


68.7370(52) 


68.01(47) 




8 


36.9335(11) 


30.17(6) 




9 


24.58756(35) 


20.87(5) 




10 


18.38542(17) 


16.294(34) 


18.437(40) 


15 


8.6540342(98) 


8.388(4) 


8.6507(57) 


20 


6.132159678(47) 


6.0645(69) 


6.1326(4) 


30 


4.29899341467(73) 


4.2919(51) 






plus sign and shrinks to zero for the minus sign when 8 — > 0(q — > 4 + ). 

Next we determine the constant a and j3 in the asymptotic form (^2|) for 
q — > 4 + . The constants a can be determined by the value of /x in the limit 
of q — > 4 + . We read from Fig. 12 that 

a = 0.073 ± 0.002 . 

On the other hand, the constant B is given by the asymptotic value of the 
following combination for n > 3: 



(24) 



Their asymptotic values should be common to each n. The behavior of the 
quantity in Eq.(p4|) for the ordered and disordered phases are given in Fig. 17 
and 18 respectively. The graphs for n = 3 and n = 4 show clear convergence 
even at 9 = giving the value of B very close to each other, that is, B ~ 0.38(3) 
We also notice that, although the results for n — 5 and n — 6 have considerably 
large errors in the small 9 region so we cannot say anything definite, the values 
of the quantity in Eq. ( p4| ) at the points where the errors begin to grow are also 
close to this value of 0.38(3). Thus we may estimate the constant B as 

B = 0.38 ±0.03 . 

In conclusion the Pade data of the large-g series for the energy cumulants con- 
vince us that the conjecture ( p2| ) for the asymptotic form is true at least for 
n = 2, 3 and 4 and suggest its correctness for n = 5 and 6. 
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Tabic 8: Estimates of the specific heat C for the ordered phase in comparison 
with the low-temperature series and the Monte Carlo simulations. 



q 


large- q 


low-temp. 


Monte Carlo 


5 


2885.8(34) 






6 


205.780(32) 






7 


68.5128(22) 


52.98(61) 




8 


36.62347(31) 


28.47(47) 




9 


24.203436(65) 


20.03(27) 




10 


17.937801(18) 


15.579(83) 


17.989(40) 


15 


7.99645871(22) 


7.721(28) 


7.999(3) 


20 


5.360768767(13) 


5.2913(53) 


5.3612(4) 


30 


3.41289525542(28) 


3.4012(10) 





Table 9: Estimates of the energy cumulants for the disordered phase (n = 
3-6). The values in the brackets are the data of the Monte Carlo simulations. 



q 


p3 
t d 


b d 


b d 


fT 


5 


2.1597(27) xlO 9 


7.4380(99) xl0 lf> 


4.90(39) xlO 22 


4.9(15) xlO 29 


6 


2.04752(87) xlO 6 


1.07367(47) xlO 11 


1.092(30) xlO 16 


1.68(38) xlO 21 


7 


9.8404(16) xlO 4 


8.2877(14) xlO 8 


1.363(14) xlO 13 


3.44(52) xlO 17 


8 


16420.4(11) 


4.69139(37) xlO 7 


2.625(11) xlO 11 


2.30(20) xlO 15 


9 


4880.80(16) 


6.69421(26) xlO 6 


1.8024(36) xlO 10 


7.67(37) xlO 13 


10 


2002.017(32) 


1.600281(33) xlO 6 


2.5155(25) xlO 9 


6.29(16) xlO 12 




[2015(26)] 


[1.583(64) xlO 6 ] 


[2.40(26) xlO 9 ] 


[5.8(18) xlO 12 ] 


15 


176.78507(17) 


32323.501(59) 


1.168612(78) xlO 7 


6.779(13) xlO 9 




[176.01(76)] 


[2.93(20) xlO 4 ] 


[1.170(59) xlO 7 ] 


[7.2(11) xlO 9 ] 


20 


54.8121419(72) 


4904.2443(16) 


870078.5(88) 


2.47889(69) xlO 8 




[54.7(4)] 


[4905(89)] 


[8.42(32) xlO 5 ] 


[2.33(16) xlO 8 ] 


30 


15.88750157(25) 


667.244592(17) 


55619.490(41) 


7.44878(13) xlO 6 



27 



Table 10: Absolute value of the estimates of the energy cumulants Fo for the 
ordered phase (n = 3-6). The values in the brackets are the data of the Monte 
Carlo simulations. 



q 




r O 


F^ 


F^ 


5 


2.1577(84) xlO y 


7.368(33) xlO 15 


5.07(21) xlO 22 


5.15(46) xlO 29 


6 


2.0467(28) xlO 6 


1.0698(22) xlO 11 


1.105(13) xlO 16 


1.775(60) xlO 21 


7 


9.8469(57) xlO 4 


8.2777(91) xlO 8 


1.3691(65) xlO 13 


3.586(52) xlO 17 


8 


16460.8(46) 


4.6939(32) xlO 7 


2.6335(57) xlO 11 


2.356(16) xlO 15 


9 


4905.23(73) 


6.7098(35) xlO 6 


1.8084(19) xlO 10 


7.799(27) xlO 13 


10 


2018.31(17) 


1.6071(11) xlO 6 


2.5268(14) xlO 9 


6.367(12) xlO 12 




[2031(26)] 


[1.55(22) xlO 6 ] 


[1.98(51) xlO 9 ] 


[2.8(14) xlO 12 ] 


15 


181.4445(19) 


32878.80(86) 


1.185486(57) xlO 7 


6.8744(12) xlO 9 




[180.67(76)] 


[2.93(20) xlO 4 ] 


[0.92(15) xlO 7 ] 


[4.5(13) xlO 9 ] 


20 


57.05595(23) 


5054.1460(86) 


892352.53(68) 


2.538205(72) xlO 8 
[1.38(26) xlO 8 ] 




[57.09(29)] 


[4.79(22) xlO 3 ] 


[7.37(70) xlO 5 ] 


30 


16.7446548(22) 


702.788233(56) 


58121.979(30) 


7.760343(17) xlO 6 
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Figure 15: The behavior of the ratio of the absolute value of the fifth energy 
cumulants \F^\ to x 11 / 2 . 
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Now, we analyze the series of the magnetization cumulants. The nth mag- 
netization cumulants M^ n > is well known to behave for (3 — > (it at q = 4 as 

MW,(-l)«M(")~Ag(0^- 2 (25) 

and we can make a conjecture parallel to the conjecture for the energy cumulants 
by Bhattacharya et ai: this relation will also be kept in the limit q — > 4 + with 
fi = fit, which implies using the asymptotic behavior of the correlation length 
in Eq.(pl|) in this limit that 

M { d n \ (-l) n M^ ~ Bfy-t n - 2 . (26) 

Therefore we use the same method of analysis as in the case of the energy 
cumulants. 

At first, we have checked that the approximants for M^C P have the best 
convergent results just &tp= (15n— 16) / 4, and next we have estimated the mag- 
netization cumulants from the Pade approximants for j^^£( 15 ™- 16 )/ 4 whose 
M and N are greater than 8 and the denominators does not possess zero at 
< z < 1/2. The results for n = 2 and 3 in the ordered and disordered phases 
are listed in Table 11 for integer values of q > 5. The estimates of the magne- 
tization cumulants have also surprisingly high accuracy. Our data are orders of 
magnitude more precise than the result of the Monte Carlo simulation by Janke 
and Kapplerfl for q = 10,15,20 (also listed in Table 11). In Figs. 19 and 20, 
we plot |M( n '|/:x;( 15 ™~ 16 ' ) / 8 versus 9 for n = 2 and 3, respectively. The constants 
B^l and B^ a in the asymptotic behavior in Eq.(|2^) might be different from 
the values at 9 = in the figures. The reason is the following. The asymptotic 
form (||) of the latent heat C for q — ► 4+ is multiplied by the correction factor 
that is an even function of #|| while the asymptotic form (||) of the magne- 
tization cumulants is plausible to be multiplied by the correction factor 
that is not an even function of 6, since in the case of the first cumulant(i.e., 
the magnetization) the correction factor is not an even function of #||. If it 
would be the case, the combination | M W y x (i5n-i6)/8 wou i(j no t b e an even 
function of 9, then it could not be approximated well arround 9 — by the 
function of z, which is an even function of 9. (On the contrary in the case of 
the energy cumulants the asymptotic form (|2^) is plausible to be multiplied 
by the correction factor that is an even function of 9, since the latent heat is 
just the difference between the two first energy cumulants. So the combination 
p{n)£in-i j g ex p ec t ec j to be well approximated by the Pade approximants in z 
even around 9 = 0.) Therefore, we will give the prediction of the values of 
based on the assumption that / £( 15 "- 16 )/ 8 varies as a linear function of 9 
around 9 = and one can obtain the values at 9 = by linear extrapolation 
from the curves around 9 ~ 0.5 in Fig. 19 and 20; 

B [ p = (2.05 ± 0.05) x 10~ 3 , B< 2) = (1.72 ± 0.02) x 10~ 3 
Bf ] = (8.0 ± 0.1) x 10~ 5 , = (8.57 ± 0.05) x 10~ 5 . 
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Table 11: Absolute value of the estimates of the magnetic cumulants M d ] Q and 
M \ ' The values in the brackets are the data of the Monte Carlo simulations. 

a,o 



q 


Mf 




5 


9.128(31) xlO 4 


9.011(23) xlO 4 


6 


628.21(37) 


666.46(39) 


7 


70.5446(94) 


77.304(13) 


8 


19.35941(71) 


21.5251(13) 


9 


8.057894(95) 


9.01050(24) 


10 


4.232764(18) 


4.738206(66) 






[4.744(42)] 


15 


0.730421382(66) 


0.80569715(69) 






[0.8052(28)] 


20 


0.3093656816(19) 


0.335564221(11) 






[0.33551(81)] 


30 


0.121320123911(18) 


0.12867774409(11) 



q 






5 


8.445(28) xlO 11 


9.336(21) xlO 11 


6 


3.3527(36) xlO 7 


3.7654(41) xlO 7 


7 


4.0328(17) xlO 5 


4.5420(35) xlO 5 


8 


2.6335(57) xlO 4 


3.3477(59) xlO 4 


9 


5116.70(50) 


5.7293(28) xlO 3 


10 


1405.906(76) 


1566.25(99) 






[1532(74)] 


15 


42.43916(31) 


46.25783(43) 






[45.3(1.0)] 


20 


7.94284(13) 


8.523886(23) 






[8.44(33)] 


30 


1.364338789(50) 


1.43676408(17) 
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0.0022 



0.002 - 




Figure 19: The behavior of the ratio the second magnetization cumulants 
to x 7 / 4 . 
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In conclusion the Pade data of the large-g series for the magnetization cumulants 
also convince us that the conjecture (|2^) for the asymptotic form is true at least 
for n = 2 and 3. 

6 Summary 

The finite lattice method was applied to generate the large-g expansion of the 
nth energy cumulants (n = 0, 1, • • • , 6) and of the nth magnetization cumulants 
(n = 1,2 and 3) at the phase transition point for the ordered and disordered 
phases in the two-dimensional (/-state Potts model. It enabled us to extend the 
series of the energy cumulants to the 23rd order in z = 1/ ^/q for n = 0, 1, 2 and 
to the 21st order for n = 3, • • • , 6 from the 10th order generated by Bhattacharya 
et al. using the standard diagrammatic method. It also enabled us to generate 
the new series of the magnetization cumulants to the 21st order. We made 
the Pade analysis of the series for the product of each cumulant and some 
powers of the latent heat that is expected to be a smooth function of z if 
the conjecture by Bhattacharya et al. is true that the relation between the 
cumulants and the correlation length in the asymptotic behavior for (3 — > (3 t at 
q = 4 will be kept in their asymptotic behavior at the first order transition point 
for q — > 4 + . It allowed us to present very precise estimates of the cumulants 
at the transition point for q > 4. The accuracy of the estimates is orders 
of magnitude higher than the results of the Monte Carlo simulations and the 
low-(and high-) temperature series for q > 7, and for instance of about 0.1 
percent at q = 5 where the correlation length is as large as a few thousand 
lattice spacings and there is no previous data to be compared with, and of a few 
percent even at the limit of q — 4 where the correlation length diverges. The 
resulting asymptotic behavior of the cumulants for q close to 4 is consistent with 
the conjecture by Bhattacharya et al. not only for the second energy cumulant 
(or the specific heat) but also for the higher energy cumulants. It is important 
that the asymptotic behavior of the magnetization cumulants is also consistent 
with the conjecture, which was not investigated by Bhattacharya et al.. This 
strongly suggests that the conjecture would be the case for all the physical 
quantities that diverge or vanish at the first order phase transition point for 
q^4+. 
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